General linewidth formula for steady-state multimode lasing in arbitrary cavities 
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A formula for the laser linewidth of arbitrary cavities in the multimode non-linear regime is derived 
from a scattering analysis of the solutions to semiclassical laser theory. The theory generalizes 
previous treatments of the effects of gain and openness described by the Petermann factor. The 
linewidth is expressed using quantities based on the non-linear scattering matrix, which can be 
computed from steady- state ab initio laser theory; unlike previous treatments, no passive cavity 
or phenomenological parameters are involved. We find that low cavity quality factor, combined 
with significant dielectric dispersion, can cause substantial deviations from the Schawlow-Townes- 
Petermann theory. 
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The intrinsic linewidth of a laser arises from quantum 
fluctuations and would be zero in the absence of spon- 
taneous emission. It is the most important property of 
lasers which arises from the quantization of the electro- 
magnetic field. Its value depends on the properties of 
the specific laser cavity and gain medium, and was first 
calculated in the seminal work of Schawlow and Townes 
(ST), who found the famous linewidth formula [1 
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where ujo is the frequency of the laser mode, j c is the 
linewidth of the relevant passive cavity resonance, and P 
is the modal output power. Note that in this formula the 
properties of the gain medium are absent. 

Improved theoretical analyses over the next several 
decades found three multiplicative corrections to the ST 
formula, all of which tended to increase the linewidth, 
in some cases by large factors [2j |3] . One correction fac- 
tor arises from incomplete inversion of the gain medium, 
and a second one from indirect phase fluctuations due 
to the instantaneous intensity change caused by spon- 
taneous emission (the Henry a factor) [3]. The third 
correction and the main focus of this Letter is the Pe- 
termann factor, K. First discovered in the context of 
transverse gain-guided semiconductor lasers [4] and sub- 
sequently generalized [2, 5tfTU], this factor arises from the 
non-Hermitian nature of the laser wave equation, due to 
the presence of the gain medium as well as the openness 
of the laser cavity (i.e. spatially non-uniform outcoupling 
loss). It always leads to an enhancement of the linewidth, 
even with uniform gain and no gain-guiding. Typically, 
it is calculated from the non-orthogonal passive cavity 
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where the integrals are taken over the cavity [IT]. In ef- 
fect, the Petermann factor changes the ST linewidth by 
the replacement 7;? — >• K^. This is a significant correc- 
tion for lasers with large outcoupling, in the range 1.1- 
1.6 for the conventional semiconductor lasers studied in 
Ref. [2 . We shall refer to the standard theory, inclu- 
sive of the Petermann factor, as the Schawlow-Townes- 
Petermann (STP) theory. 

The extensive and impressive literature on the Peter- 
mann factor [2j [Il lTOl fT2] has, with one major exception 
[9], only treated single- mode lasing near threshold, ne- 
glecting the effects of spatial hole-burning. And apart 
from a recent paper by Schomerus [12 , the literature has 
exclusively treated one-dimensional or waveguide lasers, 
and thus is not directly applicable to the wide vari- 
ety of complex laser cavities developed during the past 
twenty years, such as microdisk and deformed-disk, pho- 
tonic crystal, and random lasers. In this Letter, we de- 
rive a general formula for the intrinsic laser linewidth 
in arbitrary cavities, which is valid far from threshold, 
with strong spatial hole-burning, and in the multimode 
regime. The formula relates the linewidth to a non-linear 
self-consistent scattering matrix (5- matrix), and is based 
on the recently-developed Steady-state Ab initio Laser 
Theory (SALT) [13H16]. 

SALT is a method for solving the steady-state proper- 
ties of arbitrary lasing structures, without directly inte- 
grating the semiclassical laser equations. "Semiclassical" 
here refers to the fact that the field is treated via the clas- 
sical Maxwell equations, whereas the properties of the 
gain medium are obtained from a quantum-mechanical 
calculation of a multi-level atom. SALT treats the open- 
ness of the cavity exactly, and the non-linear modal inter- 
actions and gain saturation are included to infinite order. 
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Its results agree well with numerical integration of the 
laser equations, but it is computationally much more ef- 
ficient [17 , 18 . It has been applied to complex laser struc- 
tures such as random [15] and photonic crystal lasers [19] . 
We shall show that the quantum input-output theory of 
Refs. [20j [21] can be used to calculate quantum fluctu- 
ation properties ab initio, in terms of quantities obtain- 
able from SALT. SALT associates each laser mode with a 
scattering pole — an eigenstate of a classical nonlinear S- 
matrix with infinite eigenvalue — at a real frequency. We 
derive a formula for the linewidths of a multimode laser in 
terms of the residues of these poles and a certain norm of 
the lasing eigenstate. For a low-Q cavity, the generalized 
linewidth formula typically finds substantial deviations 
from the STP prediction: typically the linewidth is sig- 
nificantly less than the standard theory predicts, and in 
the random laser example shown below, the laser has an 
anomalous power-dependence near threshold. 
The multimode SALT equations are [16]: 
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where ^ is the fi th steady-state lasing mode, uo^ is its 
frequency, e c is the passive cavity dielectric function, 7^ 
is the gain medium linewidth, uj a is the atomic transi- 
tion frequency, D$(r) is the (possibly spatially- varying) 
pump, and T v = 7^/(7^ + (<jJ v - u a ) 2 ) is the gain 
curve. The effective pump D(f) contains an infinite- 
order nonlinear "hole-burning" term, which gives rise to 
mode competition and gain saturation in a quantitatively 
precise manner. These coupled, time-independent, non- 
linear equations are solved with the boundary condition 
of purely outgoing waves with frequency uo^ at infinity; 
the solution algorithm is discussed in Refs. [T5 ] fT6 | ITS]. 

From the solution to (|3|, we can compute a self- 
consistent ^-matrix for any complex frequency cj, not 
just the discrete lasing frequencies uj^. By definition, 
this S- matrix has one or more poles on the real-cj axis, 
at uj = uo^. It can be used to study the effects of vacuum 
fluctuations and spontaneous emission [22 . Suppose the 
cavity has scattering channels indexed by j = 1, 2, . . . , N 
(e.g. waveguide modes or spherical waves, depending on 
the scattering geometry). The input and output photon 
operators, denoted by ai, . . . , a at and 61, . . . , bjy respec- 
tively, obey an "input-output" relation [23] : 

j p 

Here the frequency Q is measured from the lasing fre- 
quency of interest, which we denote by co>o, ^ = oj — ujo. 
The dps are ladder operators for the external reser- 
voirs corresponding to the gain medium, with the index 
p denoting appropriate degrees of freedom in the cav- 
ity/reservoir [24 . 



In order for a, 6, and d to obey canonical commutation 
relations, e.g. [oj(Q), dj (£}')] = Sij 5(Q—Q f ), the ^-matrix 
must be related to the reservoir coupling coefficients by 
the fluctuation-dissipation relation [23] 

55 f - VV f = 1, (5) 

where 1 is the N x N identity matrix. Next, we define 

a j (t) = -^=Jdna j (n)e- iQt , (6) 

and similarly for bj(t) and d p (t), describing quantum am- 
plitudes for the field envelopes. Inserting into Q gives 
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The first term describes scattering of input photons, and 
the second describes emission from the gain medium. 

S is strongly constrained by its symmetries. Firstly, 
optical reciprocity [25] implies that S can be written as a 
symmetric matrix, so it has the eigenvalue decomposition 
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where each \i/j n ) denotes a right eigenvector of S with 
eigenvalue s n , and (^*| denotes its unconjugated trans- 
pose. These eigenvectors are bi-orthogonal ((VCJ^n) = 
for m 7^ n) and power- normalized ((t/ViIVVi) = !)• 

Suppose that e c (r) is real. The ^-matrix of the passive 
cavity is unitary, and for a high-Q cavity with a resonance 
near cjq, one of the eigenvalues is approximately [26] : 
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where ip is an irrelevant phase factor and j c is the cavity 
lifetime. The eigenvalue is unimodular for real £7, and, as 
required by time-reversal symmetry, its poles and zeros 
lie at conjugate positions in the complex Q plane. 

Adding gain pushes the zero and pole up in the com- 
plex frequency plane. The eigenvalue takes the form 
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where T z and T p are the distances of the zero and pole 
from the real axis. The lasing threshold is reached as 
T p — >• 0~; within the high-Q approximation the eigen- 
value takes the form (|9| with T z « 7 C (the zero moves up 
the same distance as the pole) . This approximation leads 
directly to the ST formula (high Q will imply K « 1). 
For arbitrary Q, the ^-matrix near Vt = takes the form 
([9]), with a generalized residue Yl(VL) replacing Vt — iT z . 
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We denote the eigenvector corresponding to this diverg- 
ing eigenvalue by In the S'-matrix decomposition (J8|, 
the term with sq dominates, so we can write 



so 
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Using this together with Eq. (JsJ) gives 
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This equation is satisfied by the ansatz 
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u is some vector satisfying ^2 p u p u p 
is the z th 
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where 

ty l L , is the z th component of the ^-matrix eigenvector 
for the lasing mode. Note that this relation applies not 
just to the first lasing mode at threshold, but also for 
above-threshold steady-state lasing modes, using the self- 
consistent, non-linear ^-matrix obtained from SALT. 

Inserting ([12]) into ^ and performing the resulting 
contour integrals gives 



k(t) 
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(14) 



Thus each output photon is a superposition of incoming 
photons and reservoir excitations from all earlier times. 

Above threshold, the gain medium undergoes stimu- 
lated emission, and the laser field acquires a mean value, 
B i: so that Eq. (13 becomes: 



bi{t) = Bi 
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where Bi, the steady-state classical outgoing field ampli- 
tude in channel z, is related to ty l L by 



\BA 2 = 
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where P is the total output power of the mode. 

Due to the fluctuation operator F(t), the phase of the 
output field has a quantum uncertainty; the rate at which 
this uncertainty increases with time gives the laser coher- 
ence time scale. The fluctuation-induced phase changes 
are fed back into the classical value of Bi, causing a ran- 
dom drift in the phase of the laser field. We ignore this 
feedback, instead taking a fixed value for Bi for all t. 
This is justifiable because the integrand in ( 15 ) vanishes 
exponentially for t' < — T, where T = 1/Y P will turn out 
to be the coherence time. The calculations below apply 
to times much shorter than T. 

We choose the arbitrary global phase of Bi to be real 
and positive for the specific channel z, and study the 
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FIG. 1. (color online) Output power and cavity decay rates 
7^ for two uniformly pumped one-dimensional microcavity 
lasers. A slab of gain material with background e = 1.2, 
bounded on the left by a perfect mirror and on the right by 
an e = 9 slab (5% of the total length) acting as a partially 
transmitting mirror (left schematic). A random laser consist- 
ing of 50 slabs of gain material, each with background e uni- 
formly distributed in [1, 1.2] (right schematic). Both systems 
exhibit two-mode lasing at the high end of the pump range. 
Plots (a),(b) show modal output powers vs. the normalized 
pump 16 . Plots (c),(d) show the square of the generalized 
cavity decay rate ^\ = |r£| 2 /|J/J^i,| 2 which determines the 
linewidths according to Eq. pi] ). Solid and dashed curves 
denote modes 1 and 2 respectively. The horizontal dotted 
lines show the conventional result, if 7c, computed from the 
passive cavity quasimodes, which fails for the random laser. 



quantum fluctuations of the phase via the Hermitian 
quadrature operator [28 
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which for small phase angles corresponds to the phase of 
the laser output in channel z. Using (14) and (15), we 
compute the quantity (^(^1)^(^2))? taking (a) = \d) = 
(a\(ti) CLjfe)) = and taking the white noise correlator 



(d p {t 1 )di{t 2 ))=f p 5 pv 5{t 1 -t 2 ), 



(18) 



where f p = [P2/CP2 — Pi)]p describes the local popu- 
lation inversion [23]. The zero-point contributions to 
(Oi(ti)0i(t2)) from the photon input and the gain medium 
cancel exactly, leaving 
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where P is the modal output power from Eq. (H6j), and 
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FIG. 2. (color online) Laser linewidth vs. inverse modal out- 
put power 1/P, for the two lasers studied in Fig. [I] The 
linewidths are computed using Eq. (21), assuming the inver- 
sion factor / = 1. (a) The high-Q cavity laser linewidths 
show the standard 1/P dependence for both modes, (b) 



The linewidth of the first mode of the random laser devi- 
ates strongly from the 1/P Schawlow-Townes-Petermann de- 
pendence at lower pump values. At large pump values the 
linewidths of both mode 1 (solid curve) and mode 2 (dashed 
curve) vary as 1/P, but with values roughly half that of the 
standard STP prediction (dotted curve). 



is the inversion factor correction mentioned at the begin- 
ning of this Letter. 

The phase uncertainty accumulated over time At is 
([Oi(t + At) - 6i(t)} 2 ) = AujAt + 0(At 2 ), where 
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This is our central result: a general linewidth formula in 
which \T L \ 2 /I^^lI 2 = j1 replaces the quantity Kj 2 in 
the conventional Schawlow-Townes-Petermann linewidth 
formula. We can think of as a generalized cavity decay 
rate, corrected for the presence of gain, openness, hole- 
burning, and saturation. It is calculable ah initio, with no 
phenomenological parameters, from the non-linear classi- 
cal ^-matrix of SALT. The lasing eigenvector is found by 
diagonalizing the S'-matrix at each lasing pole, and the 
residue Y^ is found by numerically integrating the rele- 
vant eigenvalue of the ^-matrix around the pole. Eq. (21 ) 



only includes the contribution to the laser linewidth from 
direct phase fluctuations; the indirect phase fluctuations 
[3] have been omitted for simplicity. 

The relation of the Petermann factor to the residue 
of the lasing pole for a waveguide laser was emphasized 
early on by Henry [6 , and developed for more general 
cavities in an ^-matrix formulation in Refs. [12, 29 , but 
in all previous cases for a single lasing mode at thresh- 
old, i.e. without non-linear effects. Goldberg et al [9] 
gave an excellent and detailed analysis of the linewidth 
for multimode lasing, including non-linear effects, but us- 
ing an approach applicable only to one-dimensional cav- 
ities with spatially uniform dielectric functions. To our 
knowledge, our Eq. (21), combined with SALT, is unique 



in providing a quantitative method for calculating the 
intrinsic laser linewidth in arbitrary cavities and pump 



profiles in the multimode, non-linear regime. Assuming 
steady-state multimode lasing exists, the present theory 
makes no significant further approximations, and hence it 
can be used to evaluate the validity of the STP linewidth 
formula 121. 



We can connect Eq. (21) to previous results involv- 



ing quasi- modes, such as Refs. [U [5j [7], by examining 
the ^-matrix of a passive cavity. A quasimode ip(r) is 
a purely-outgoing solution to the wave equation for a 
passive cavity with dielectric function e(r), at complex 
frequency uj P1 where Im(co> p ) = — 7 c /2. Let \I/ be the S- 
matrix eigenvector for this pole, normalized by ^Sfr = 1, 
and let Y be the residue of the eigenvalue. It can be 
shown that 
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Here the spatial integrals are taken over the cavity. For 
real e(r), and in the limit \k p \ ^> Y ~ j c , (22)-(23) give 
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Thus, in this slightly generalized version, Kj 2 is approx- 
imately equal to our 7 ?, w hen evaluated for the passive 
cavity. Note that both ( 24 ) and its active-cavity general- 
ization in Ref. [12] do not include the effects of dielectric 
dispersion, which can have a significant effect on 7^. 

Fig. [I] compares 7^ to Kj 2 for two one-dimensional 
microcavity multimode lasers: a high-Q, uniform cav- 
ity for which the two quantities agree rather well, and 
a low-Q random laser, for which major deviations are 
found. For the random laser, at pump strengths up to 
four times threshold, 7^ for the first lasing mode depends 
strongly on P, causing the overall power dependence to 
depart substantially from the standard 1/P dependence 
(Fig. [2|. For higher pump strengths, 7^ is approximately 
constant, but the conventional linewidth prefactor Kj 2 
overestimates it by almost a factor of two. In the stan- 
dard theory, the STP linewidth is expected to be a lower 
bound set by field quantization, but insofar as the usual 
STP formula relies on passive cavity quantities it is not 
a reliable bound. Analysis of our results indicates that 
this deviation from the STP theory arises from low cav- 
ity Q and from the frequency dispersion of the dielectric 
constant of the gain medium, which cause a significant 
reduction of the residue of the lasing pole at threshold 
compared to its value in the passive cavity. We do not 
believe that the apparent violation of the STP bound in- 
dicates any new quantum fluctuation properties of the 
laser. In future work, our generalized linewidth formula 
will allow such issues to be studied systematically. 
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